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Recent work in light flavour hadron physics is reviewed. In particular, I discuss the signifi- 
cance of the progress achieved with light dynamical quarks on the lattice for the effective low 
energy theory of QCD. Also, I draw attention to some puzzling results from NA48 and KTeV 
concerning the scalar form factor relevant for A" M 3 decay - taken at face value, these indicate 
I ' physics beyond the Standard Model. 
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J> ' 1 Introduction 

OO 

At low energies, the most important characteristic of QCD is that the energy gap is very small 
(f) ■ - a consequence of the fact that the u- and d-quarks happen to be very light. In the theoretical 

\£2 ■ limit where the quark masses m u and rrid are set equal to zero, the Hamiltonian of QCD acquires 

an exact, spontaneously broken symmetry. The spectrum of the theory does then not have an 
energy gap at all: in the limit, the pions become massless particles, playing the role of the 
Goldstone bosons which necessarily occur if a continuous symmetry spontaneously breaks down. 

The remarkable theoretical progress made in light flavour hadron physics in recent years 
relies on the fact that these properties can be used to construct an effective theory ("chiral 
perturbation theory", referred to as XPT), that allows us to analyze the low energy structure 
of QCD in a controlled manner. In this framework, all of the Green functions formed with the 
quark currents can be calculated in terms of the coupling constants occurring in the effective 
Lagrangian, order by order in the chiral expansion. For recent reviews of XPT, I refer to 
Scherer , Bijnens 2 and Colangelo. 3 An up-to-date account of our knowledge of the effective 
coupling constants can be found in a recent conference report by Ecker. 4 

In my talk, I focused on a few selected aspects of this development. In particular, I empha- 
sized the fact that the progress achieved in the numerical simulation of QCD on a lattice made 
it possible to reach sufficiently light quarks, so that the extrapolation to the quark masses of 
physical interest can be done in a controlled manner, using XPT. I will discuss this in some de- 
tail below. Another recent development which I will briefly report on, concerns the semileptonic 
decay K — > 7t/xia Both KTeV and NA48 have recently published new results on the scalar form 
factor of this decay, which are in flat conflict with a venerable low energy theorem, established 
by Callan and Treiman in the sixties. 5 If these results are confirmed, then the Standard Model 
is in conflict with observation in one of those reactions which we thought are best understood. 

For lack of space, I cannot cover the third topic which I dealt with at these Rencontres: the 
progress made in understanding the properties of the interaction among the pions. This inter- 



action plays a crucial role in many contexts - the Standard Model prediction for the magnetic 
moment of the muon is perhaps the most prominent example. To close this introductory section, 

1 briefly list the corresponding keywords and indicate where more information about this can 
be found. The irir scattering amplitude has been calculated to NNL in XPT. 6 The resulting 
representation is very accurate in the interior of the Mandelstam triangle, but in the physical 
region, the convergence of the series is rather slow. The range of energies where the chiral 
representation yields meaningful results can be extended with the inverse amplitude trick, but 
one is then leaving the territory where model independent statements can be made. There is a 
general method that does not suffer from such shortcomings: in a limited region of the complex 
plane, dispersion theory imposes a set of exact relations between the real and imaginary parts 
of the partial wave amplitudes, the Roy equations. 7,8 The region where these equations are valid 
includes the poles on the second sheet generated by the lowest resonances of QCD: a, p, u, 
/o(980), ao (980). The crucial parameters that control the low energy properties of the scatter- 
ing amplitude in this framework are the two subtraction constants. It is convenient to identify 
these with the two S-wave scattering lengths, because the low energy theorems of XPT make 
very sharp predictions for these. Together with the low energy theorems, the Roy equations pin 
the scattering amplitude down within remarkably small uncertainties. 9 The angular momentum 
barrier ensures that the S- and P-waves dominate at low energies, but the framework also yields 
very accurate predictions for partial waves with higher angular momenta. On the basis of this 
method, it can be demonstrated beyond reasonable doubt that the lowest resonance of QCD 
carries the quantum numbers of the vacuum and the position of the corresponding pole can be 
worked out rather accurately. 8 For an overview of these developments, I refer to the proceedings 
of Chiral Dynamics 2006 and the references quoted therein. 10 

2 Size of the energy gap 

In order to illustrate the progress achieved on the lattice, I consider one of the key issues in 
QCD: understanding the size of the energy gap. The quark masses m u and are very small, 
but they are different from zero. Accordingly, the Hamiltonian of QCD is not invariant under 
chiral rotations - the quark mass term breaks chiral symmetry and there is an energy gap: 
the symmetry breaking equips the Goldstone bosons with a mass. The quark masses m u , rrid 
represent a quantitative measure of the strength of the symmetry breaking. We know that the 
symmetry breaking is very small, but, unfortunately, the Standard Model does not offer an 
understanding of why that is so - the entire fermion mass pattern looks bizarre and yet remains 
to be understood. 

As pointed out by Gell-Mann, Oakes and Renner, 13 the square of the pion mass is propor- 
tional to the strength of the symmetry breaking, oc (m u + m^). This property can now be 
checked on the lattice, where - in principle - the quark masses can be varied at will. Fig. 1 shows 
the result for two recent lattice simulations of QCD with two flavours. In view of the fact that 
in these calculations, the quarks are treated dynamically, the quality of the data is impressive. 
The masses are sufficiently light for XPT to allow a meaningful extrapolation to the quark mass 
values of physical interest. The results indicate that the ratio M%/(m u + md) is nearly constant 
out to values of m u ,m^ that are about an order of magnitude larger than in nature. 

3 Lattice determinations of the effective coupling constants £3 and £4 

The Gell-Mann-Oakes-Renner relation represents the leading term in the expansion in powers 
of the quark masses. At next-to-leading order, this expansion contains a logarithm: 





(1) 




M. Liischer, Lattice 2005 11 Ph. Boucaud et al. [ETM Collaboration] 

Figure 1: Lattice results for M% as a function of the quark mass 



Chiral symmetry fixes the coefficient of the logarithm in terms of the pion decay constant F, 
but does not determine the scale A3 of the logarithm. A crude estimate was obtained more than 
20 years ago, 14 on the basis of the SU(3) mass formulae for the pseudoscalar octet: 

A 2 

0.18 GeV < A 3 < 2GeV £ 3 = In = 2.9 ± 2.4. (2) 

The logarithmic term implies that the lines in Fig. 1 cannot be straight. For the central value, 
£3 = 2.9, the formula (1) yields only little curvature, but if £3 was at the lower (upper) end of 
the quoted range, the plot of Ml versus m u = would be strongly bent upwards (downwards), 
visibly departing from the lattice results shown in Fig. 1, already at the lowest quark mass value. 
Evidently, these results strongly constrain the value of the coupling constant £3 . 



Table 1: Determinations of the effective coupling constants £3 and £4 





XPT 9 ' 14 


MILC 17 


Del Debbio et al. 18 


ETM 12 


£3 


2.9 ±2.4 


0.6 ± 1.2 


3.0 ±0.5 


3.62 ±0.12 


£4 


4.4 ±0.2 


3.9 ±0.5 




4.52 ±0.06 



The first row in table 1 compares the number in equation (2) with recent determinations 
of £3 on the lattice. The second row lists the analogous results for the coupling constant £4, 
which controls the quark mass dependence of the pion decay constant. To my knowledge, the 
first lattice calculation of effective coupling constants based on dynamical quarks was carried 
out by the MILC collaboration. 15 In that project, the coupling constants L4, L5, Lq, Lg, which 
occur in the effective chiral SU(3)xSU(3) Lagrangian at first nonleading order, were determined 
by analyzing the quark mass dependence of M n , Mr, F n and Fk by means of XPT. The corre- 
sponding values of £3 and £4 are readily worked out, using standard one loop formulae. 16 The 
numbers quoted in the third column are obtained from a recent update of the MILC analysis, 
which relies on staggered quarks. 17 The results obtained by Del Debbio et al. 18 are based on two 
flavours of Wilson quarks, while the European Twisted Mass Collaboration 12 uses two flavours 
of twisted mass Wilson quarks. All of the numbers are consistent with the values found in XPT, 
but some of them are more accurate. Unfortunately, the results obtained with staggered quarks 



do not agree well with those based on Wilson quarks. Possibly, this is related to the fact that 
the latter calculations treat all quarks except u and d as infinitely heavy - as emphasized by 
Stern and collaborators, the strange quark may play a significant role here. 19 Lattice results for 
3 flavours of light Wilson quarks should become available soon, so that it should be possible to 
identify the origin of the difference. 




The values of the low energy constants £3 and £4 also enter the theoretical prediction for the 
7T7T scattering lengths. This is illustrated in Fig. 2, where theory is compared with experiment 
(the lattice results for £3, £4 are converted into corresponding values for a(j,ao using XPT and 
dispersion theory 9 ' 10 ). While the lattice results, DIRAC and the NA48 data on the cusp in 
K — > 3tt confirm our predictions for a°)> a 0> the K e ^ data of NA48/2 give rise to a puzzle: 
the phase extracted from the transition amplitude deviates from the theoretical prediction for 
<5g — 5{ . The discrepancy originates in the fact that neutral kaons may first decay into a pair of 
neutral pions, which then undergoes scattering and winds up as a charged pair. As pointed out 
by Colangelo, Gasser and Rusetsky, the mass difference between the charged and neutral pions 
affects this process in a pronounced manner: it pushes the phase of the transition amplitude up 
by about half a degree - an isospin breaking effect, due almost exclusively to the electromagnetic 
interaction. The dash-dotted line in Fig. 2, which is taken from the talk of B. Bloch-Devaux at 
KAON 2007, shows the likelihood contour (x 2 = Xmin + 2.3) of the so corrected, preliminary 
NA48 /2 data. The intersection with the region allowed by the low energy theorem for the scalar 
radius yields a{j = 0.220(9). However, there is a discrepancy with the E865 data, for which 
the likelihood contour is shown as a dashed line (kindly provided by G. Colangelo): the isospin 
correction spoils the good agreement between these data and the prediction. The discrepancy 
only concerns the region M nn > 350 MeV. While E865 collects all events in this region in a 
single bin, the resolution of NA48/2 is better. The fit to all K e ^ data is therefore dominated by 
NA48/2. For a detailed discussion of these issues, I refer to the talks by B. Bloch-Devaux, G. 
Colangelo and J. Gasser at KAON 2007. I conclude that the puzzle is gone: K e ^ confirms the 
theory to remarkable precision. 



4 Puzzling results in K^s decay 



The low energy theorem of Callan and Treiman 5 predicts the size of the scalar form factor of 
the decay K — > TTfiv at one particular value of the momentum transfer, namely t = M\ — M 2 : 

f (M 2 K - Ml) = + 0(m u ,m d ) . (3) 

Within QCD, the relation becomes exact if the quark masses m u and rrid are set equal to zero. 
The corrections of first nonleading order, which have been evaluated long ago, 20 are tiny: they 
lower the right hand side by 3.5 x 10~ 3 . In the meantime, the chiral perturbation series of fo(i) 
has been worked out to NNL. 21 As pointed out by Jamin, Oiler and Pich, 22 the curvature of 
the form factor can be calculated with dispersion theory, so that the prediction for the value at 
t = M\_ — M 2 can be converted into a rather accurate prediction for the slope: Ao = 0.0157(10). 
The dispersive representations of Jamin et al. 22 and Bernard et al. 23 agree very well: theory 
reliably determines the curvature of the form factor. 

Very recently, the NA48 collaboration published their analysis of the K® s form factors. 2 
Their result for the scalar slope, Ao = 0.0117(7)(1), is in flat conflict with the prediction of 
Jamin, Oiler and Pich. Using the parametrization proposed by Bernard et al., 23 NA48 obtains 
/o(M£ - M 2 ) = 1.155(8)(13) x /+(0). The value of F K /F n is sensitive to the theoretical input 
used for /+ (0) . Accurately measured branching ratios 25 and the value of V u d (which has been 
determined to high precision on the basis of nuclear (3 decays) imply Fx/F n = 1.241(4) x /+(0). 
So, if the NA48 results are correct, then the second term on the right hand side of (3) must 
amount to a contribution of 0.086(17) x /+(0). For a correction of 0(m u , m^), this size is unheard 
of (as mentioned above, the one loop approximation for this term amounts to -0.0035 and is thus 
smaller by a factor of about 20). At the current experimental accuracy, the radiative corrections 
may play a significant role. For decay, these are known to one loop of XPT 26 but the data 
yet remain to be analyzed in this framework. 

The NA48 experiment is not the first to measure the slope of the scalar form factor. The first 
results were obtained in the seventies. In particular, the high statistics experiment of Donaldson 
et al. had confirmed the theoretical expectations with a slope of Ao = 0.019(4). More recent 
experiments, however, came up with quite different results. In particular, three years ago, the 
KTeV collaboration at Fermilab arrived at a remarkably small scalar slope: Ao = 0.01372(131). 
Analyzing 0.54 million charged kaon decays, ISTRA, on the other hand, obtained a much higher 
value: A = 0.0 196(12)(6). 29 If the results of NA48 as well as those of ISTRA were correct, then 
the dependence of the form factors on the momentum transfer would have to show very strong 
isospin breaking - that would be extremely interesting in itself. 

My conclusion is that the experimental situation calls for clarification. In particular, an 
analysis of the charged kaon decays collected by NA48 might help removing the dust. There 
are not many places where the Standard Model fails. Hints at such failures deserve particular 
attention." 
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"In the meantime, a preliminary analysis of the KLOE data on the scalar K® 3 form factor became available. 25 
The result for the slope, Ao = 0.0156(26), is in excellent agreement with the theoretical prediction, thus removing 
the puzzle ... at the price of generating a new one: KLOE is in conflict with NA48. 



References 



1. S. Scherer, Eur. Phys. J. A 28 (2006) 59, hep-ph/0512291. 

2. J. Bijnens, Prog. Part. Nucl. Phys. 58 (2007) 521, hep-ph/0604043. 

3. G. Colangelo, Lectures given at the Int. School of Subnuclear Physics, Erice, Italy (2006). 

4. G. Ecker, in Effective Theories of Colors and Flavors from EURODAPHNE to EURIDICE, 
Kazimierz, Poland (2006), to be published in Acta Phys. Polon. hep-ph/0702263. 

5. C. G. Callan and S. B. Treiman, Phys. Rev. Lett. 16 (1966) 153. 

6. J. Bijnens, G. Colangelo, G. Ecker, J. Gasser and M. E. Sainio, Phys. Lett. B 374 (1996) 
210; Nucl. Phys. B 508 (1997) 263; ibid. B 517 (1998) 639 (E). 

7. S. M. Roy, Phys. Lett. B 36 (1971) 353. 

8. I. Caprini, G. Colangelo and H. Leutwyler, Phys. Rev. Lett. 96 (2006) 132001. 

9. G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B 603 (2001) 125. 

10. H. Leutwyler, in Proc. 5 th International Workshop on Chiral Dynamics, Theory & Exper- 
iment, Durham/Chapel Hill, USA (2006), hep-ph/0612112. 

11. M. Liischer, PoS LAT2005 (2006) 002, hep-lat/0509152. 

12. Ph. Boucaud et al. [ETM Collaboration], hep-lat/0701012. 

13. M. Gell-Mann, R. J. Oakes and B. Renner, Phys. Rev. 175 (1968) 2195. 

14. J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

15. C. Aubin et al, [MILC Collaboration], Phys. Rev. D 70 (2004) 114501. 

16. J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465, section 11. 

17. C. Bernard et al., in Proc. 5 th International Workshop on Chiral Dynamics, Theory & 
Experiment, Durham/Chapel Hill, USA (2006), hep-lat/0611024. 

18. L. Del Debbio et al, JHEP 0702 (2007) 056 and JHEP 0702 (2007) 082. 

19. S. Descotes-Genon, N. H. Fuchs, L. Girlanda and J. Stern, Eur. Phys. J. C 34 (2004) 201. 

20. J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 517. 

21. P. Post and K. Schilcher, Eur. Phys. J. C 25 (2002) 427. 
J. Bijnens and P. Talavera, Nucl. Phys. B 669 (2003) 341. 

V. Cirigliano, G. Ecker, M. Eidemuller, R. Kaiser, A. Pich and J. Portoles, JHEP 0504 
(2005) 006. 

22. M. Jamin, J. A. Oiler and A. Pich, JHEP 0402 (2004) 047, hep-ph/0401080. 

23. V. Bernard et al, Phys. Lett. B 638 (2006) 480. 

24. A. Lai et al. [NA48 Collaboration], hep-ex/0703002. 

25. M. Palutan, talk presented at KAON 2007, transparencies accessible via LNF 
Frascati/Conferences/21-05-07 Kaon International Conference. 

26. V. Cirigliano, M. Knecht, H. Neufeld, H. Rupertsberger and P. Talavera, Eur. Phys. J. C 
23 (2002) 121. 

27. G. Donaldson et al, Phys. Rev. D 9 (1974) 2960. 

28. T. Alexopoulos et al. [KTeV Collaboration], Phys. Rev. D 70 (2004) 092007. 

29. O. P. Yushchenko et al., Phys. Lett. B 581 (2004) 31 (value converted to units of M n +). 



